Computation of Grobner basis for systematic 
encoding of generalized quasi-cyclic codes 

Vo Tarn Van Hajime Matsui Seiichi Mita 
Dept. Electronics and Information Science 
Toyota Technological Institute 
2-12-1 Hisakata, Tenpaku, Nagoya, 468-8511, Japan 

Abstract 

Generalized quasi-cyclic (GQC) codes form a wide and useful class 
of linear codes that includes thoroughly quasi-cyclic codes, finite ge- 
ometry (FG) low density parity check (LDPC) codes, and Hermitian 
codes. Although it is known that the systematic encoding of GQC 
codes is equivalent to the division algorithm in the theory of Grobner 
basis of modules, there has been no algorithm that computes Grobner 
basis for all types of GQC codes. In this paper, we propose two al- 
gorithms to compute Grobner basis for GQC codes from their parity 
check matrices: echelon canonical form algorithm and transpose algo- 
rithm. Both algorithms require sufficiently small number of finite-field 
operations with the order of the third power of code-length. Each al- 
gorithm has its own characteristic; the first algorithm is composed of 
elementary methods, and the second algorithm is based on a novel 
formula and is faster than the first one for high-rate codes. Moreover, 
we show that a serial-in serial-out encoder architecture for FG LDPC 
codes is composed of linear feedback shift registers with the size of the 
linear order of code-length; to encode a binary codeword of length n, 
it takes less than 2n adder and 2n memory elements. 

Keywords: automorphism group, Buchberger's algorithm, division 
algorithm, circulant matrix, finite geometry low density parity check 
(LDPC) codes. 



1 Introduction 



Low density parity check (LDPC) codes were first discovered by Gallager [9] 
in 1962 and have recently been rediscovered and generalized by MacKay [18] 
in 1999. The methods of constructing LDPC codes can be divided into two 
classes: random construction and algebraic one. 

Random constructions of irregular LDPC codes [5J [18] [22] have shown 
the performance near to the Shannon limit for long code lengths of more 
than 10 7 bits. On the encoding of random LDPC codes, Richardson et al. 
[23] proposed an efficient encoding method with the decomposition of the 
generator matrix into low triangular matrices, which was improved by Kaji 
[TTj and Maehata et al. [20] with another triangular (or LU-) factorization. 
Both methods of encoding are based on the matrix multiplication. 

There are many algebraic constructions of LDPC codes [lj[6j[7j[12j[26j, 
which belong to a class of quasi-cyclic (QC) codes and provide efficient de- 
coding performance. Another remarkable algebraic construction of LDPC 
codes is finite geometry (FG) codes [32] [15]; These codes are divided into 
Euclidean (or affine) geometry (EG) codes, which are included in QC codes, 
and projective geometry (PG) codes, which are included not in QC codes 
but in broader generalized quasi-cyclic (GQC) codes (cf. Figured]). It can 
be stated briefly that GQC codes increase the randomness for QC codes and 
vary each length of cyclic parts in QC codes. 

For several classes of QC LDPC codes, Fujita et al. [8] proposed efficient 
encoding with circulant matrices and division technique. With regard to 
GQC codes, which includes the algebraic LDPC codes, Heegard et al. [TO] 
showed that the systematic encoding was equivalent to the division algorithm 
of Grobner bases, which generalize the generator polynomials in cyclic codes. 
According to this work, Chen et al. [3] constructed an encoder architecture. 
Thus, the encoding problem for GQC codes was changed into the computa- 
tion of Grobner basis. For the computation of Grobner basis for encoding 
GQC codes, Little [16] provided an algorithm for Hermitian codes, and Lally 
et al. [33] provided an algorithm for QC codes. However, there has been no 
algorithm applicable to all GQC codes. 

In this paper, we propose two algorithms for computing the Grobner 
bases, which encode GQC codes, from their parity check matrices. The 
first algorithm is based on Gaussian elimination, and the second algorithm 
is the generalization of Lally et aVs algorithm. Both algorithms employ 
Buchberger's algorithm to create a Grobner basis from codewords. Moreover, 
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Figure 1: Inclusion-exclusion relation of various linear codes. 



in order to show its efficiency, we prove that the number of circuit elements 
in the encoder architecture is proportional to code-length for finite geometry 
codes. 

A part of the first proposed algorithm to compute Grobner basis was 
already known to some specialists in coding theory. Kamiya et al. [12] 
announced that an encoder was obtained with fundamental row operation for 
a QC LDPC code from Euclidean geometry. Recently, Little [T7| announced 
a similar result for a Hermitian code. Our object is to provide algorithms 
computing the Grobner bases for all GQC codes even in the case requiring 
column permutation. On the other hand, the second proposed algorithm is 
based on a novel formula that produces Grobner basis from that of the dual 
code. The special case of our formula was found by Lally et al. [14J for 
QC codes. In order to extend it to the case of GQC codes, we provide our 
formula with a completely different proof from that of Lally et al.'s formula. 
Both algorithms have 0(n 3 ) order of the computational complexity, where 
n is the code-length, and in fact for high-rate codes, we can show that the 
second has less complexity than the first. 

Although the size of the encoder architecture for general GQC codes 
exceeds the linear order of code-length because of the number of orbits (cyclic 
parts), Chen et al. [3] proved that it had the linear order for Hermitian 
codes. We newly prove that it also has the linear order for FG codes. While 
Richardson et al.'s and Kaji's methods for general LDPC codes run by the 
linear order of finite-field operations, our encoder architecture for FG codes 
can achieve not only the linear order of operations but also the linear order of 
circuit elements and no latency. In addition, our encoder architecture for the 
binary FG LDPC codes requires only adder elements without multiplication 
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(i.e., no AND element). 

This paper deals with all GQC codes; Siap et al. |24J mainly focused 
on one-generator GQC codes. Another example of GQC codes is the class 
of algebraic geometry (AG) codes with automorphism groups [10], including 
Hermitian codes [16]. It is worthy to notice that GQC codes include two 
remarkable classes of Hermitian codes and some PG codes outside QC codes. 
Thus, GQC codes form the vastest algebraic class in linear codes that holds 
compact encoder architecture. Therefore, we can choose more appropriate 
and high-performance codes from GQC codes than those from QC codes. 

This paper is organized as follows. Section 2 provides the definition of 
GQC codes and the techniques of Grobner basis. Section 3 provides the 
details of the first: echelon canonical form algorithm. Section 4 provides 
the details of the second: transpose algorithm. In section 5, we estimate 
the computational complexity of proposed two algorithms. In section 6, we 
prove the linearity of the circuit-scale of the encoder architecture for FG 
LDPC codes. Finally, we conclude this paper in section 7. 

2 Preliminaries 

Throughout the paper, we denote A := B if A is defined as B. First, we de- 
scribe the definition and module structure of generalized quasi-cyclic codes. 
Then, we review Grobner basis of modules over polynomial rings; the com- 
plete theory of Grobner basis is referred to [2] [4] , and that of automorphism 
group and orbit is referred to [19J. 

2.1 Definitions 

Consider a linear code C C F™ of length n, where q is a prime power and 
F g is g-element finite field. Let S be the set of locations (that is, coordinate 
positions) of codewords in C: C 3 c = (c s ) s€ s- Without loss of generality, we 
set S = {1, 2, ■ • • , n}. Now suppose that there is a decomposition of S, 

m m 

S = \JO t , \S\ = n = J2 li > l i : =\°i\> (!) 
i=i i=i 

and accordingly, decompose any codeword c G C into m shortened codes: 

C = (Ci, C2, • • ■ , c m ), (2) 
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where q is a shortened codeword dropping components outside . Consider 
simultaneous local cyclic shift a of each q satisfying 

o-(c) := (cr(ci), ■ ■ • ,cr(c m )), 

a(ci) := (cj^_i, c i)0 , • ■ ■ , e iA _ 2 ) (3) 
for a = (ci t o,Ci t x, ■ ■ ■ ,c it i.-i). 

Definition 1 If we have m < n and er(c) G C for all c G C, then we call a 
pair of C and o a generalized quasi-cyclic code (GQC code). □ 

If C is GQC, thus we obtain a nontrivial a of the automorphism group 
Aut(C) of C. Conversely, if Aut(C) includes a ^ 1, then the cyclic group (a) 
generated by a defines orbit O(s) := {o~ l (s) \ a 1 G (a)} of s G S. Note that 
we have O(s) = O(s') for s' G O(s), and that we have 0(s) = {s} if a(s) = s. 
Then, S is equal to the disjoint union of distinct orbits as described in ([[]), 
where O t := O(sj) for some Sj G S 1 . We can regard a as the simultaneous 
shift (E]) of {q}. Thus, we have shown that the class of GQC codes agrees 
with the class of linear codes with nontrivial Aut(C) D (a). 

Remark 1 Thus, we see that each d := {c^} decides a cyclic code. How- 
ever, in general the whole C does not agree with the combined code Yl Ci 3 
[ci • • • c m ], since the individual shift, for example (a(ci), c 2 , ■ • ■ , c m ), does not 
generally belong to C. We will see the difference between C and f| Cj at (j7|) 
in subsection 12.31 

Remark 2 Siap et al. [24J define GQC codes as Fjtj-submodules of a certain 
module M at (jl]) in the next subsection, where we will see that their definition 
of GQC codes is equivalent to our definition. 

Usually, the generator matrix of a linear code indicates the matrix whose 
rows are linearly independent and compose a basis of the linear space. We 
often relax this definition for convenience; we call a generator matrix of a 
linear code the matrix whose rows are not too many and contain the basis. 

Example 1 Consider the linear code C\ C Fj defined by a generator matrix 
as below. 
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Since the second row plus the third row equals the fourth row, we see that 
the dimension of C\ is three. If we apply the permutation a given by ([3D, 
then a codeword in C\ is transferred into another codeword in C\. Thus, C\ is 



made from 4 cyclic codes defined by (1), (111 
(and two all-zero codes), and C\ is a GQC code with 3 orbits 
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Note that, if l\ — h — ■ • • — l m , then C is a quasi-cyclic code [T4"j [T5] [2T] . 
Moreover, if m = 1, then we come back to cyclic code. In order to increase 
the randomness of the codes, it is desirable that we can combine various 
circulant matrices (cf. Section @J to generate new GQC codes. The code in 
Figure [2], a) shows us a 2-orbit GQC code constructed from four matrices. 
On the other hand, Figure [21 b) shows that the code is also obtained from 
four matrices but is not a GQC code. 



a) 



b) 



Figure 2: Intuitive models of generator matrices made from four matrices. 
Model a) defines a 2-orbit GQC code, but Model b) does not define a GQC 
code. 



2.2 Module structure of generalized quasi-cyclic codes 

Let C be a GQC code with a permutation a. Under the action of (a), we 
can decompose c G C into m shortened codes as described in (j2J). Pick c; and 
label it as Cj = (cij) where j = 0, ■ ■ • , U — 1 with U := \Oi\. For convenience, 
we decide that the second index is an integer modulo Zj, and the permutation 
cr of (j3J) means c(cjj) = Ci t (j-i m odk) for allz = 1, ■ • • , m and j — 0, • ■ • , Zj — 1. 
Then, a codeword in C can be represented as an m-tuple of polynomials in 
¥ q [t): 

C = (Ci(t),C 2 (t), ■ • • ,Cm(t)), 
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k-l 

where q(t) = c i,j^ ■ Thus, C is regarded as a linear subspace of M, where 

j=0 

m 

M := (F 9 M/ (** - 1) ) (4) 

t=i 

and F g [t]/ (i^ — l) is the quotient ring by an ideal (t h — l) := (t li — l) ¥ q [i\. 
Moreover, we can regard the action of a as the multiplication of t as follows: 

\% — 1 \% 1 \% 1 

tci(t) = ^ ti+1 = r -i |/; = a{r >j)' J - ( 5 ) 

j=0 j=0 j=0 

where "=" means the equality modulo (t li — l) . We can see that multiplying 
c by t is equivalent to permuting the codeword locally cyclically by a. Thus, 
C is closed under the multiplication by t and C is considered as an ¥ q [t\- 
submodule of M. For convenience to compute Grobner basis, we consider 
the following natural map: it : ¥g[t] m — > M. Let be the i-th standard basis 
vector in Fjtj-module F g [t] m , that is, 

e x := (l,0,0,---,0), e 2 :=(0,l,0,---,0), 

•■■ , e m := (0,0,0,--. ,1), 

and Xi := (t lt — l) e« for i = 1, • ■ • , m. Define C := 7r _1 (C), which is a 
submodule of F g [t] m and is generated by all codewords in C (regarded as 
vectors in F g [t] m ) and all X^s, that is, 

C = C+ (X, I i = I,-- ■ ,m) , (6) 

where (Xi \ % — 1, • • ■ , m) indicates the submodule generated by all X;. 

2.3 Grobner basis of ¥ q [t] -module 

We call an element of the form Pei a monomial in F g [t] m . Then, any polyno- 
mial vector in F g [t] m can be represented as a linear combination of monomials. 

Although Grobner basis of a submodule in F g [t] m is determined for each 
monomial ordering, only the following two orderings are required in this 
paper. The position over term (POT) ordering [10] on F 9 [t] m is defined by 

>pot t k ej if i < j, or i = j and I > k. Then, we have e\ >pot &2 >pot 
■ • • >pot e m . 
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Similarly, the reverse POT (rPOT) ordering is defined by t l ei > r poT t k ej 
if i > j, or i = j and I > k. Then, we have e\ < r poT ^2 <rPOT • • • <rPOT e m . 

For a polynomial fit) := Yli=o fit* e with fd 7^ 0, we define degree 
of f(t) by deg(/) := d, and we say that / is monic if fd — 1. Thus, we can 
define two types of Grobner bases for the submodule C of F g [i] m associated 
with an m-orbit GQC code C. 

Definition 2 We define POT Grobner basis of C as the following set Q = 
{gi, 92, ■ ■ ■ , 9m} of polynomial vectors 

9i = (9u(t), 9i2(t), ••• , 9i m (t)), 

» = (o, te(t ), , *.(*)), (7) 

9m = (0, • ■ • , 0, g mm {t)) 

such that gi, ■ ■ ■ , g m £ C and gait) has the minimum degree among the vectors 
of the form (0, • • • , 0, Cj(i), • • • , c m {t)) G C with 7^ 0. If ga 's are monic 
and Q satisfies deg < deggjj for all I < i < j < m, then we call it reduced 
POT Grobner basis. Moreover, we define rPOT Grobner basis of C as the 
following set Ti. = {hi, h 2 , ■ ■ ■ , h m } of polynomial vectors 

h 1 = (h u (t), 0, ••• , 0), 

: : : (8) 

it), , h m —\ m —i 

(t), 0), 

hm \hml\t), j h mm —\it),h mm it)) 

such that h\, ■ ■ ■ , h m 6 C and hu(t) has the minimum degree among the 
vectors of the form (ci(t), • ■ ■ , Cj(t), 0, • • • ,0) EC with Cj(t) 7^ 0. // ha 's are 
monic and Ti, satisfies deg h%j < deg hjj for all 1 < j < i < m, then we call 
it reduced rPOT Grobner basis. □ 

Since (t li - l) e^'s are included in C, the diagonal polynomials ga and ha 
divide r* — 1, and both Grobner bases of C exist. If = for all i 7^ j, then 
C agrees with the combined code Y[ C% as noticed at Remark [TJ 

From any Grobner basis, we can easily obtain the reduced Grobner basis 
by fundamental row operations of polynomial matrix. Each GQC code has 
its unique reduced Grobner basis. 
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The Grobner basis of C has two important roles; one is that it generates C, 
and the other is division algorithm (which is stated in the next subsection). 
Any element c G C has the following expression 

c = P 1 (t)g 1 + --- + P m (t)g m , (9) 

where Pi(t) G F 9 [i]. If deg(P;(fe) < k for all i, then we have c e C strictly. 
For the use of encoding, we define redundant monomial as Vei with < j < 
deg ga(t) (standard monomial in [10J). The other types of monomial Pei 
with deg ga(t) < j < k are called non-redundant (or information) monomial. 
It follows from (P) that the number of information monomials equals the 
dimension of C. 

By minimizing deg ga, we can obtain Q from the generator matrix. This 
procedure is called Buchberger's algorithm, which is now described for our 
situation. 

Buchberger's algorithm 

Input: A k x n generator matrix G of a GQC code C 
Output: A POT Grobner basis Q = {gi, ■ ■ ■ ,g m }. 
Step 1. Regard the row Gi of G as polynomial vector Gi 
W q [t] m . 

Step 2. For j = 1 to m; 

gj ■= J2i=j QijGi with {Qjj, ■ ■ ■ , Q kj } such that 

9jj '■= gcd{Gjj, ■■■ ,G kj } = J2i =j QijGij . 
For I — j + 1 to k; 

Gi := Gi - Ri(t)gj with R t := Gij/gjj. 

Step 3. If gj is zero vector, then put gj := Xj. Reduce Q by fundamental 
row operations. □ 

Example 2 Consider again Example HJ Buchberger's algorithm is applied to 
G\. We obtain g\ = (1, 1 + t, 1) by adding the third row to the first row. 
Since the other polynomial vectors are the multiple of gi, we obtain gi = X2 
and #3 = A 3 as shown in Figure [31 Then, the information monomials are 
t 2 ei, tei, e\ and the redundant monomials are t 2 e 2 , te 2 , e 2 , e 3 . □ 
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gl = (1, 1 + t, 1), 

g 2 = (0, 1 + t 3 , 0), 
g 3 = (0, 0, 1 + t). 
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Figure 3: The reduced POT Grobner basis of C\ and monomials. 



2.4 Systematic encoding algorithm 

Once a Grobner basis Q = {gi, • • • , g m } of C is obtained, then division algo- 
rithm with respect to Q can be applied to u G F g [t] m to obtain the following 
representation 

u = Qi (t)gi H h Q m {t)g m + u, (10) 

where Qi(i) G ¥ q [t], and u = with degMj < deg^j. In 

other words, u is a unique linear combination of redundant monomials. It 
follows from ([9]) and (|T0l) that m G C <^ U = (0, • • • , 0), which generalizes the 
condition of codewords in cyclic codes. Then, the encoding of C is described 
as follows. 

Systematic encoding algorithm 

Input: Information symbols «6Fj and Grobner basis Q = {gi, ■ ■ ■ , g m }. 
Output: Encoded codeword c G C. 

Step 1. Calculate u G Fjt] m as a linear combination of information symbols 

and information monomials. 

Step 2. Put u x = (u n (t), ■ ■ , u lm (t j) := u; 

For i — 1 torn; 

Find Qi(t) and Ui(t) such that 

uu{t) = Qi(t)g u (t) + Ui(t), deg Hi < deggu . 
Calculate u, + i := Ui — Qi(t)gi (g M) 

= (Ux(t),--- ,Ui{t),Ui + x ii+l {t), ■ ■ ■ ,Wi+l, m (t)). 
Put « := («i(t), ■ ■ • = «- Eili^i^)^ in M - 

Step 3. By subtraction c := u — u, we obtain the encoded codeword c G C. 
□ 
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Step 2 itself is called division algorithm, which generalizes the classical 
polynomial division in the encoding of cyclic codes. Thus, another merit of 
considering the reduced Grobner basis is that it reduces the computational 
complexity of the division algorithm. 

Example 3 We reuse the GQC code C\ with the reduced Grobner basis Q\. 
The information symbols for C\ can be taken as the coefficients of information 
monomials {ei, te±, t 2 ei}. We apply the systematic encoding algorithm to 
encode message (1, 0, 1). First, we put u := t\ + • tt\ + t 2 e\ = (1 + t 2 , 0, 0). 
Then, we divide u by Q\ to obtain the remainder (or parity symbols) u: 

u = u- {l + t 2 ) gi = (0,t + t 2 ,0), 

where the last equality follows from (0, 1 +t 3 , 0) = (0, 0, 0) in M. Since u con- 
tains only redundant monomials, we finish the division algorithm. Thus, the 
encoded polynomial vector is u — u = (1 + 1 2 , t + 1 2 , 0) and the corresponding 
encoded codeword c is equal to (1010110). Since we have 



(1 1 1 1 0) = (0 1 1 0) 
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we can check that c G Ci, which is required. □ 

3 Computing Grobner basis from parity 
check matrix with echelon canonical form 

In this section, we consider the problem about computing Grobner basis, 
which generates a GQC code, from a given parity check matrix. In many 
situations, each GQC code C is specified by a parity check matrix, that is, 
the generator matrix of its dual code C . Since we have Aut(C) = Aut(C _L ), 
both codes are viewed as the submodules of the same M at (p!J. 

Before describing the proposed algorithm, we remind that elementary row 
operations can be used to simplify a matrix and we obtain echelon canonical 
form [2TJ, which is defined as follows: 

• Every leftmost non-zero value is 1. 
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Figure 4: Outline of computing Grobner basis for encoding by echelon canon- 
ical form algorithm from parity check matrix. 



• Every column containing the leftmost non-zero value has all zero other 
entries. 

• The leftmost non-zero value in any row is on the right of that in every 
preceding row. 

For example, we consider the following two matrices: 



/10100001 \ 
110 10 
10 1 

o o o o no o 
\oooooon/ 



/ 1 1 1 1 1 o o \ 
110 10 
10 10 1110 
110 10 10 
\ 1 1 1 1 1 ) 



The left matrix is the echelon canonical form of the right matrix. 

By using the echelon canonical form, we can compute the Grobner basis 
of C from parity check matrix. The flow of our first algorithm is presented 
in Figure @] and described as follows: 

Echelon canonical form algorithm 

Input: Parity check matrix if of a GQC code C. 

Output: POT Grobner basis Q of C. 

Step 1. Transform H to echelon canonical form H' by Gaussian elimination. 
Step 2. Select permutation r satisfying H\ := t(H') = [I\A], and then 
G X >.= [-A T \1\. 

Step 3. Compute generator matrix G = r" 1 (G ! i). 

Step 4. Obtain Q by Buchberger's algorithm from G. □ 
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In the step 2 of the above algorithm, G\ satisfies the equation G\ x = 
and r _1 (G?i) x {r- 1 (H 1 )} T = r -1 (G?i) x H' T = 0. Therefore, we permute 
the column vectors of G\ by r _1 to obtain a generator matrix G of the GQC 
code C. 

Example 4 Let C2 be a GQC code defined by the following parity check ma- 
trix H2: 



H 2 := 
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We can see that C2 has locally cyclic property with the column permutation 
a = (1 • • • 6)(7 • • • 12) (13 • • • 15), where, e.g., (1---6) indicates permutation 
1— > 2 — > 6 — > 1, and C2 has 3 orbits: l\ = l 2 = 6 and l 3 = 3. Firstly, we 
use Gaussian elimination to transform H2 to the equivalent echelon canonical 
form 
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(11) 



If we choose the column permutation r such that the set of column location 
(1, 2, • • • , 14, 15) is mapped by r to 

(1,2, 7,8,9,10, 3,4,5,6, 11,12,13,14,15), 

then matrix ffTTl) is transformed to the standard form matrix [I\A]: 
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(Note that, in this case, r has no relation to the orbit decomposition.) Then, 
we permute the corresponding matrix [— y4 T |i] by r _1 to obtain the generator 
matrix G2 of Ci- 



( 
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By using Buchberger's algorithm, we can compute the reduced POT Grobner 
basis {91,92,93} of generator matrix G 2 : 



9i 




' 1, 


0, 1 


92 




0, 


1 + t + t 3 + t 4 , 1+t 


93 




0, 


0, 1+t + t 2 



□ 

Although this example is binary, our algorithm can be applied to all parity 
check matrix H of F^-entries. We consider in section the computational 
complexity of our algorithm to obtain Q from H. 



4 Transpose formula for POT Grobner basis 

In this section, we propose another algorithm to compute the Grobner basis 
from parity check matrix. This novel algorithm uses transpose formula (j20|) 
that is given at Theorem [2j Although Theorem 1 is not necessary for our 
computation except a scalar product (Tl3|) and Corollary [TJ we describe it for 
completeness; Theorem 1 provides the orthogonal property with respect to 
the scalar product for arbitrary Grobner bases of GQC codes. 

Firstly, we define a circulant I x I matrix as a square I x I matrix such 
that each row is constructed from the previous row by a single right cyclic 
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shift. Then, we can represent the circulant I x I matrix 



.4 



/ CLq CLi 

Q-i-l a o 



ai-i \ 
ai-2 



as a polynomial a(t) = + ait + • ■ ■ + a^it 1 ' 1 in module F g [i]/ (t l — l). 

Proposition 1 Let a(t) and b(t) represent the corresponding polynomials of 
circulant matrices A and B of size I x I, respectively. 

(i) Transpose of A is a circulant matrix corresponding to polynomial a(t) := 
ao + ai-it + ■ • • + ait 1 " 1 in module Fjt]/ (t l — l) . 

(ii) Matrix product AB equals BA and corresponds to polynomial a(t)b(t). 
In particular, we have AB = if and only if a(t)b(t) = mod (t l — l). 



(Hi) If [Ai--- A m ] * 



Bl 



holds, where Ai and Bi are circulant I x I 



matrices, then we have the corresponding polynomial J2 a i(t)bi(t) = mod 

i=i 

(f-i). 

(iv) The product 



A 



A 



* \B T ■ ■ ■ £? T ] equals a circulant matrix of size ml x ml 



m—l 



an 



d corresponds to polynomial a(t)b(t) t d mod (t ml — l) . 



8=0 



Proof: Proposition l.(i)-(iii) are easy to prove and we refer to [14j. Propo- 
sition l.(iv) can be proved by executing matrix multiplication 



.4 
A 



[ B T ■ ■ ■ B T ] 



AB T 
AB T 



AB T 
AB T 



:i2) 



From (i) and (ii), we see that the matrix (|T2l) is a circulant matrix of size ml x 
ml, and moreover, the circulant matrix AB T can correspond to polynomial 
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a(t)b{t) in module ¥ q [t}/ (t l - l). Therefore, matrix flTJ) can be represented 
by the following polynomial 

m— 1 

a(t)b(t) (1 + t l + • ■ • + = a(t)b(t) tU 

in module ¥ q [t}/ (t ml - l). □ 

Next, we define a scalar product of polynomial vectors u — (ui, ■ ■ ■ , u m ), 
v = (vi, ■ ■ ■ ,v m ) e M as 

(u,v) := J^(t)%(t) ^ ^ mod (t l - 1) , (13) 

i=l k=0 

where / is the least common multiple (1cm) of Zj's that correspond to non-zero 
Ui and Vi. We denote [u] as the matrix representation of polynomial vector 
u = (ui, ■ ■ ■ ,u m ) G M by shifting locally cyclically I times. Since U divides 
/ and t li is regarded as 1, we can represent matrix [u] by non-zero circulant 
matrices [ttj] and zero matrices. For example, assume u = (1+t, 0, 1+t 2 ), v = 
(1 + t + t 3 , 1 + t, 1 + t) e M, where h = 5, l 2 = 4, l 3 = 3, q = 2. Since 
I = lcm(Zi, l 3 ) = 15, the matrix representation of polynomial vector u agrees 
with 



"110 0" 




" 








" 




"10 1" 


110 


















1 1 


110 


















1 1 


1 1 


















1 1 


1 1 


















1 1 


110 


















1 1 


110 


















1 1 


110 


















1 1 


1 1 


















1 1 


1 1 


















1 1 


110 


















1 1 


110 


















1 1 


110 


















1 1 


1 1 


















1 1 


_ 1 1 


















_ 1 1 _ 
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We see that matrix [u] can be decomposed into two non-zero circulant ma- 
trices 

/ 1 1 \ 
110 
110 and 
1 1 
\ 1 1 / 

which correspond to polynomials (1 + t) mod (t 5 — 1), and (1+t 2 ) mod 
(t 3 — 1), respectively. The scalar product (u, v) agrees with 




(i + t)(i + 1 2 + 1 4 ) t5k + (i + 1 2 ) (i + 1 2 ) Yl t3k 



k=0 



k=0 



= (l+t 2 + t 8 + t 10 + t u + t 14 ) mod (t 15 - 1) . 

With these preparations, we can obtain the following orthogonality between 
a Grobner basis of a GQC code and that of its dual. 

Theorem 1 Let Q = {gx, • • • , g m } and TC = {hi, • • • , h m } be a Grobner basis 
of a GQC code C and that of C with respect to any ordering, respectively. 
Then, we have (gi, hj) = (hj,gi) = mod (t l — l) for all 1 < i,j < m. 

Proof: (hj,gi) = follows from (gi,hj) = easily. Shifting the component 
9i = {9a, ■■■ ,9im) locally cyclically / times, we obtain j^ -1 ^.} 
that correspond to I polynomial vectors as follows: 



9n 


9i2 


9im 


tga 


%2 


t9im 




1 9i2 


■■ t l ~ x a- 



\ 



(14) 



Let [g^ denote the matrix corresponding to I vectors (|14|) . We can represent 
[gi] by circulant matrices [gij] as follows: 



[9i 



( [9a] U 
\ [9u] U 



[9im] ^ 
[9im] J 



(15) 
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/ l9ik] \ 

where \ is the I x If. matrix made from non-zero matrix [gik] or only 
\ [ftfc] / 

from zeros. Since t 5 gi G C for all S, every rows of [g^] are codewords in C. 
Similarly, the corresponding matrix representation of Ttj agrees with 



[hi) 



( [hjl] M [hjm] ^ 

\ [hji] [h j2 ] [hjm] J 



and every rows of [hj] are codewords in C^. The relation c* (c^) T = 0, where 



c G C and cr 1 G C x , corresponds to [hj] 
have the following equivalent equation 



for all Therefore, we 



/ [9a] \ 




( [M ^ 


T ( U \ 


( M \ 












\ hi] J 






V U J 


\ M / 




f [9im] \ 







+ ••• + 



0. 



(16) 



\ [dim] J \ [hjm] ) 



If Qik = or hik = 0, then the fc-th term of (1161) is zero matrix of size I x I. 
Otherwise, by Proposition HJ(iv), the fc-th term of ffTBI) is a circulant I x I 
matrix and the corresponding polynomial agrees with 

g ik (t)h jk (t) (t l - 1) / (fr - 1) mod (t l - 1) . 

Therefore, the corresponding polynomial of ffTBI) is obtained as follows: 

i/h-i i/h-i 
gu(t)hji{t) tSh + 9a{t)h j2 {t) t&h 



5=0 



<5=0 



l/lm-1 



+ ■■■ + g im {t)hj m {t) t&lm = (9i, hj)=0 



5=0 



modulo (t l — l), which leads the theorem. 
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□ 



By using Theorem [TJ we can compute the Grobner basis Q from 7i. 
However, the computation is not straightforward because of the ambigu- 
ity "mod (t l — l)." Little et al. [16] obtained strict equalities for POT and 
rPOT diagonal components (gi, hi), which we applied to finite geometry codes 
in [27]. Now, we remove all modulo conditions. For later use, we derive a 
corollary from the argument at f[T5"j) . 

Corollary 1 Let 7i = {hi, h 2 , ■ ■ ■ , h m } be a Grobner basis ofC L , and u G M 
a polynomial vector. Then, it holds that (hi,u) = for all 1 < i < m if and 
only if u corresponds to a codeword in C. □ 

In the case of cyclic codes, if we know the generator polynomial h(t) of 
the dual code C and a(t)h{t) = t n — 1, then that of C is the reciprocal 
polynomial t dega a(i -1 ) of a{t), which agrees with t dega a(t) mod (t n — 1). 
We generalize this relation to GQC codes. Assume that Tl is rPOT Grobner 
basis of an m-orbit GQC code C^. Since Ti is a basis of C L (as described at 



Q), there exists m x m polynomial matrix A 



satisfying 



.4 







' t h - 1 





h 2 




t h - 1 













h m 




o 


t lm - 1 



It is easy to observe that A = 
(hij), namely, = if i < j. 
A(hij) = (hij) A as noticed in [1 



is a lower triangular matrix similar to 
If C 1 - (or C) is a QC code, then we have 
, but in general not commutative. We can 



calculate a™ recursively as follows: 



t k - 1 



-1 



if j > i, 
if j = i, 



M S=j+1 



(17) 



— ^ a is h sj if j < i- 



It is important fact that, if TC is the reduced rPOT Grobner basis, then 
A = (a,ij) has the similar property, that is, dega^ < degau for all i > j. Now 
we prove this by induction on j. The first step dega^-i < dega^ follows 
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from a^j-i/ii-i^-i + a^/i^-i = from (fT7j) . Suppose induction hypothesis 
deg < deg for j + 1 < 8 < i. From (IT7j) . we obtain 

deg a^- = deg ^ Oij/i^ J - deg hjj 

< max {deg a^ + deg h$j — deg /ij-,} 

j<S<i 

< deg a iiy 

which proves the fact. 

From now on, we assume that 7i is the reduced rPOT Grobner basis. We 
define transpose polynomial matrix of A by 



an 


021 


"ml 




' h ' 





0-22 













Q"mm 




b 



(18) 



where G"ij IS calculated in ¥ q [t}/ (t h - l), and not in ¥ q [t}/ (ft - l). bmce a,-,- 
is the j-th component of a polynomial row vector, it might seem natural to 
calculate in F q [t]/ {t lj — l). Nevertheless, we consider a^ modulo {t lt — l), 
which is justified by deg a^ < deg an < U and is a characteristic of GQC codes 
that is disappeared in the case of QC codes. 

The latter half of the following theorem provides the objective formula of 
POT Grobner basis for GQC codes. 

Theorem 2 Polynomial vectors (|T8|) satisfy 

Moreover, Q = {g\, ■ ■ ■ ,g m }, where 

9ij ■= r mod (t h - 1) , (20) 

determines a POT Grobner basis of GQC code C (usually not reduced). 

This formula fl20|) generalizes that of cyclic codes and that of QC codes 
by Lally-Fitzpatrick [T3j to the case of GQC codes. In [13], their formula is 
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proved by the fundamental row operation of a polynomial matrix; this proof 
cannot be applied to our case because of the complication to different orbit 
lengths. We first show ffT9l) directly from (fTTj) . then we conclude by degree 
argument. 

Proof of Theorem^ From the definition (fT3l) . it obviously holds that 
(hi,bj) = for i < j, and then we concentrate on the case of j < i. Con- 
sider two polynomial vectors hi = (ha,-- - ,ha,0,--- ,0) 6 W and bj = 
((),••-, 0,a 



hi 



1 Q>mj ) , 



where 1 < j 1 < i < m. From ffTTj) and a — a, it 
is trivial that (hk, bk) = t lk — 1 for all 1 < k < m. Thus, we may prove only 
(hi, bj) = for all j < i by induction on i—j. We denote /3y := lcm(Z.,, • • ■ , Zj), 
then (hi,bj) is computed as follows: 



(hi, = (t^ - 1) £ 



hikQ>kj 

t l *-l' 



If z — j = 1 , we have 



(^> = (t^-l) 
_ 1 



j<k<i 



th - i t'i - i 



30 



13 J 



If % — j — 2, we have 

6,-) = (t^ - 1) 

= - 1) 

From f|T7|) . we obtain 



hijCLjj hij-^-yClj+ij hiiOj- 

F 7 r ~ 

h 
h 



1 ' t^+i - 1 ' £'* 
hi t j-\-\&j-\-\j dij 



(I , I 



Therefore, (hi,bj) corresponds to 



hjj Qj%% 



{o>ijhjj + a^j+i/ij+ij + CLuhij) — 0. 



(21) 
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Suppose induction hypothesis {he, bg) = for all 9 < 8 with 5 — 9 < i — j . 
From (I17p . we receive, for all j + 1 < k < i — 1, the following equations 



^ i— l ^ fc 

/iifc = — '^^O'ishsk, a-kj = t — cbkehgj. 

5=k n H 9=j+l 



Consider the following partial summation of fT2T|) : 

i-i , i-i ^ /i-i \ ( k 

aiih H Y = Y t^~\ Y aishsk Y akehe i 

k=j+l k=j+l \S=k / \e=j+i 

= ^0 J i5hsk a k9hgj 

j+l<k<i-l k<6<i-l 

j+i<e<k 



a^/iflj ^ -hsk^ke- (22) 



j+l<0<<5<i-l 6Kfc<<5 

For all j + 1 < 6* < <5 < 2 — 1, we have {hg,bs) = by induction hypothesis. 
Therefore, from (ITTj) . the double summation (1221) is equal to 

i— 1 , i—1 

E llggagg _ 

a i9' l 8j7f g ^ ~~ 2 -> CL i9' l Sj ~ ~~ a ij' l jj ~~ a ii tl ij- 

6=3+1 9=j+l 

From this result, {hi, bj) is equal to 



- 1) 



hijCljj Q>ij h^j hfaCLij 



t l i — 1 an hjj t lt — 1 







Therefore, we have {hi,bj) = for all 1 < i ^ j < m. 

The rest of the proof is to show that (1201) determines a POT Grobner 
basis in the meaning of Definition [21 From Corollary P, we have bi G C, then 
gi G C. Thus, we may prove only that gu has the minimum degree among 
the vectors of the form (0, • • • , 0, Cj(t), • • • , c m {t)) EC with Ci(t) ^ 0. We 
notice that Ci := {q | (0, • ■ • , 0, q, • • ■ , c m ) G C} defines a cyclic code. Since 
the generator polynomial of the dual code is ha, that of Ci is gu, then gu 
has the minimum degree. □ 

By Theorem [21 we obtain the second algorithm for computing Grobner 
basis Q of m-orbit GQC code C from the parity check matrix as follows. The 
flow of this algorithm is presented in Figure 
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Buchberger's Inductive Transpose and 



Parity check 


algorithm 


Reduced rPOT 


formula 


Polynomial 


cyclic shift 


Grobner basis Q 


matrix H 




Grobner basis H 




matrix A 




generating Q 



Figure 5: Outline of computing Grobner basis by transpose algorithm from 
parity check matrix. 



Transpose algorithm 

Input: Parity check matrix if of a GQC code C. 
Output: POT Grobner basis Q of C. 

Step 1. Compute the reduced rPOT Grobner basis 7i by Buchberger's 

algorithm from matrix H. 

Step 2. Calculate A = (a^) by ffTTjl . 

Step 3. Obtain Q = {gi, ■ ■ ■ ,g m }, where 



9i — {9ij)l<j<m, 9ij 





t dcga "a 



if % > j, 
,, if i < j- 



□ 



Remark 3 We can construct the generator matrix G of GQC code C from 
its reduced POT Grobner basis {gi, ■ ■ ■ ,g m } as follows: 



G 



9u 


912 


9lm 


tgu 


tg 12 


tgim 


t xi gn 


t Xl gi2 ■■■ 


t xi gim 





922 


92m 





tg 22 


t92m 





t X2 g 2 2 ■■■ 


t X2 92m 








9mm 








t9mm 








f%m n 

ymm 


j 9im, 


and Xi := I 


i - deg g u 



\ 



(23) 



where g t = (0, • • • ,0,^. 
Since the diagonal components ga all lie in different position, the rows of 
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this matrix are linearly independent. Moreover, the total number of rows 

m 

equals ^2 (k ~ deggu) = k. Therefore, the matrix (|23|) provides the generator 
i=i 

matrix of GQC code C, which generalizes the representation for quasi-cyclic 
codes in HU. □ 



Example 5 We demonstrate the transpose algorithm. 
GQC code with l\ 



Let C3 be a binary 



= 6, l 2 




6, h 




4 defined by 






/o 1 


1 


1 


1 


10 10 





\ 


1 


1 


1 


1 


10 10 








1 1 





1 1 





10 1 








1 


1 


1 


1 


1 


1 





1 1 





1 1 





1 





10 


1 


1 


1 


1 


10 





10 


1 


1 


1 


1 


10 





10 


V 1 1 





1 1 





1 





1 / 



We calculate the reduced rPOT Grobner basis TC3 
C% by Buchberger's algorithm: 



{hi, ti2, hs} of dual code 



" hi ' 




1 + 1 6 , 


0, 





h 2 




t + t 2 + t 4 + t 5 , 


l+t 2 , 





h 




l+t + t 3 + t 4 , 


1, 


1 



There exists a polynomial matrix A 
we can calculate A inductively: 



"VI. 



satisfying A[hi] = 0. From ifPTj) . 



.4 



1, 

t + t 2 + t 3 , 
l+t 2 , 



0, 

l+t 2 + t\ 
l+t 2 , 






l + t 4 



The transpose polynomial matrix of A turns into 



an, a 2 i, a 3 i 

0, 022, Ct32 

0, 0, a 33 



1. 

0. 

0. 



t 3 + t 4 + t 5 , 
l + t 2 + t\ 
0, 



l+t 2 
l + t 2 
l + t 4 



According to Theorem [21 a POT Grobner basis of GQC code C3 can be 
computed by := t degau ciji mod (t li — l). After reduction, we obtain the 
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reduced POT Grobner basis Q 3 = {gi, g 2 , gz}'- 



' 9i ' 




' 1, 


1 + t + t 2 , 


t + t 3 


92 




o, 


l+t 2 + t\ 


1 + t 2 


93 




o, 


o, 


1 + t 4 



To check the correctness of (plj) . we calculate the generator matrix G 3 of C 3 

by m>- 



G 3 = 



/ 1 


1 


1 


1 











1 1 \ 


1 





1 


1 


1 








10 10 


1 








1 


1 


1 





10 1 


1 











1 


1 


1 


10 10 


1 


1 











1 


1 


10 1 


1 


1 


1 











1 


10 10 





1 





1 





1 





10 10 


V o 





1 





1 





1 


10 1/ 



We observe that G 3 x Hj = 0, as required. □ 

Remark 4 It should be noted that Theorem [T] is valid not only for POT and 
rPOT ordering but also for any ordering. We demonstrate Theorem [1] to the 
term over position (TOP) ordering [10] on Fq[t] m defined by t l ti >top t k ej if 
/ > k, or / = k and i < j. The reduced TOP Grobner basis {g[, g' 2 , g' 3 } of C 3 
turns into 







1 + t, 


t + t 2 + t 3 , 


1 + t 2 






1, 


1 + t + t 2 , 


t + t 3 


9' 3 




1 + t + t 2 , 


o, 






It is easy to check that (g[, hj) = for all □ 



5 Estimation of algorithms 

In this section, we estimate the computational complexity of two algorithms 
and compare the one with the other. We represent the numbers of additions, 
subtractions, multiplications, and divisions in ¥ q as the coefficients of k, A, fi, 
and v, respectively. We aim for an asymptotic estimation; we denote / ~ g 
if f/g tends to 1 as the variable tends to oo. Once we obtain / ~ g, then 
it follows that the usual notation / = O(g), which means / < eg for some 
constant c > 0. 
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First, we describe it with respect to the Gaussian elimination. Although 
it is well-known that the complexity is 0(n 3 ), we calculate the order up to 
constant factor. We can assume that a given (n — k) xn parity check matrix 
H is transformed into [I\A] with permutation r = 1 since the permutation 
costs no finite-field operation. Without loss of generality, we can assume that 
the (1, 1) component of H is non-zero. Then, dividing the other component 
of the first row by this value takes (n — l)u. Moreover, subtracting the 
multiple of the (i, 1) component and the first row for 2 < % < n — k takes 
(n — k — l)(n — l)(X+fj,). Summing up these manipulations for n — k columns, 
we obtain 

fc+i 

{iv + k)i(X + /i) + (n - i)k(X + //)} 

i=n—l 

where the last term (n — i)k(X + /i) comes from the back substitution. We 
ignore the first term iu since it contributes square order |(n + k)(n — k — 1) 
of n. Then, we obtain (u + X + /i) times 

-(n — k — l)(n — k — —)(n — k) + k(n — k) 2 , 

which is asymptotically \{n — k) 3 + k(n — k) 2 . 

Next, we describe the computational complexity with respect to the Buch- 
berger's algorithm for a given k x n generator matrix G to obtain a POT 
Grobner basis. The estimation is similar to the above; now the algorithm 
is based on polynomial gcd computation. Without loss of generality, we 
can assume that the (l,h) component of G is not zero. Then, dividing the 
other component of the first row by this value, and moreover, subtracting 
the multiple of the (i, l±) component and the first row for 2 < % < k takes 
(n — l)v+(k — l)(n — l)(A+/i). The second stage of these manipulations takes 
(n — 1)v + k(n — 2) (A + /i) since the first row has a polynomial whose degree 
is greater than that of the other rows. Summing up these manipulations for 
the first orbit, we obtain 



h-i 

E 

j=di 



j + J2k) {u + k(X + fi)} 



i=2 



-(n-l)(A + /x), 



where we denote di := degtfa and the last term, which is ignored because of 
less contribution, comes from the special situation stated above at the first 
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row. Furthermore, simplifying this and summing up for all orbits, we obtain 
(v + X + fx) times 

ll±^+f; fcV-j + l), (25) 

i=j+l J 

which is bounded by nk 2 , since the second bracket < n. 

It is necessary that we estimate the reducing computation of Grobner 
basis, which corresponds to the back substitution of polynomial matrix. For 
POT Grobner basis, we must start from reducing g 12 and not gi m . The length 
of <7i2 is I2 in the worst case, and we have to eliminate I2 — d% values. Thus, 
the reduction of gn takes (I2 — (I2) (g?2 + Y4L3 h) (A + //). Summing up for all 
Qij (i < j), we obtain (A + ji) times 

m / m \ 

5> - \ d i- E MO'" 1 )' (26) 
i=2 \ i=i+i / 

which is bounded by rank. 

For the total complexity of Buchberger's algorithm to obtain the reduced 
Grobner basis, we must add (l2"fiT) to (1231) . Since we can bound (12"6"|) by the 
summation of j from 1, the last bracket of fl25|) is changed into /c. Then, we 
observe that the total complexity is still nk 2 . 

The final stage of estimation is to calculate the number of operations 
required for computing the polynomial matrix A from the polynomial matrix 
(hij) by ffT7|) . It should be noted that the multiplication of two polynomial 
a(t) and b(t) requires (degadeg5)ft + (1 + dega)(l + deg&)/z operations, and 
that the division of a(t) by b(t) requires deg 6(deg a — deg b)(\+/i) operations. 
We denote ej := deg ha; then we have 

m m 

deg an = h- €i, 22 e * = k ' E ^ ~ e d = n ~ k - 

i=l i=l 

From f|T7|) . we see that the computation of Oy (j < i) is separated into 
two steps: Yl a ishsj an d its division by hjj. Since dega^ < U — e« and 
deg /i^j < €j, the complexity of the first step is bounded by 



E ~ e i~ !)( e i _ !) K + (*< _ e i) e jA i 

5=j+l 



(27) 
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Since we can check by direct calculation that the coefficient of k in (|27l) is 
bounded by (i—j)(li—€i)ej, thus ( 1271) is bounded by (i— j)(h— ei)ej(K+fi). The 
second step requires — ei)ej(X+fi). Hence, the complexity of computing 
is bounded by (i — j + l)(Zj — ej)ej(/c + A + //). On the other hand, we see from 
(ITT)) that the complexity of computing is bounded by {U — e^e^A + //), 
which is viewed as the case of i = j for (z — j + — £i)ej(A + //). Summing 
up these results, we obtain (re + A + //) times 

- j + l)(7i - ei)ej < mk(n - k). 

i=l j=l 

Therefore, the complexity of computing polynomial matrix A is estimated as 
mk(n — k). 

Thus, we obtain estimation formulae. 

Echelon canonical 1. ,. 9 , 9 

: -(n - kY + k(n - kY + nk 2 
form algorithm 3 

Transpose algorithm : n(n — k) 2 + mnk + mk(n — k) 

We can observe that both algorithms have the same rough order 0(n 3 ) of 
computational complexity. 

For the comparison of two algorithms, we assume m = |n, |n, |n to elim- 
inate m in the estimation formula for transpose algorithm. Since there exist 
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special GQC codes that satisfy m = n — 1, these assumptions are not valid 
for general GQC codes. However, for effective GQC codes, these assumptions 
are reasonable; actually, FG LDPC codes have further less number of orbits 
than these assumptions (cf. Table [1] in the next section). Figure [6] is the 
comparison between the coefficients of n 3 in the above estimation formulae 
under the assumptions, where the curves near 1 represent limits. Thus, we 
can conclude that, for the effective high-rate GQC codes, the computational 
complexity of the transpose algorithm is far less than that of the echelon 
canonical form algorithm. 

6 Estimation of the circuit 

In the previous sections, we have proposed two algorithms to calculate the 
reduced Grobner basis of the form (JTj) that generates m-orbit GQC code C. 
In [3], Chen et al. have developed a serial-in serial-out hardware architecture 
to encode information symbols systematically with POT Grobner basis as 
an application of results in Heegard et al. [HI]. The architecture generalizes 
classical encoder of cyclic codes and consists of division circuits by gait) and 
multiplication circuits with gij(t) {i < j). 

We quote the estimation of their hardware complexity from [3]. The 
total numbers of finite-field adder elements A m and memory elements (shift 
registers) D m are given as follows: 

m m—1 m 

A m < ^deggu + ( de S^i + 1) 

i=l i=l j=i+l 

m—1 

< (n — k) + '^^(m — i) deggu < m{n — k), 
i=i 

m m—1 m m—1 

D m < ^deggu + 2J de S#u + /S^ + X ) - m ( n ~ k) + k, 

i=l i=l j=i+l i=l 

where 5i := max (ki — 2, k 2 — 2, ■ - • , ki_ x — 2,ki — 1), and ki := U — deggu. 
We can conclude that the hardware complexity for GQC codes is nearly 
proportional to the code length since m is small compared to n. 

For more practical estimation, we focus on the finite geometry (FG) 
LDPC codes [13j[15j[21j as an important class of GQC codes. There are 
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two types of FG LDPC codes: type-I and type-II. Type-I FG LDPC codes 
are defined by the parity check matrix composed of incidence vectors (as 
rows) of lines and points in finite geometries (Euclidean geometry (EG) and 
projective geometry (PG)) and are cyclic codes. Type-II FG LDPC codes 
are defined by the transposed parity check matrix of type-I and are not cyclic 
but GQC codes. Therefore, we concentrate on type-II FG LDPC codes. We 
quote the required properties of this type of codes from [25J. We denote n' 
and k' as the corresponding values of type-I codes. 

1. li < l 2 — ■ ■ ■ — l m (Actually, it becomes the equality for EG codes.) 

2. gn = ■■ ■ = g m -i,m-i = 1 and degg mm = n-k 

3. ( n -k) = (n' - k') < n' = l m 

The last two properties follow easily from the fact that the dual codes of 
FG LDPC codes are the one-generator GQC codes, which corresponds to the 
case of I — 1 in [25], eq.(23)]. Therefore, the reduced POT Grobner basis 
Q = {gx, ■ ■ ■ , g m } of type-II FG LDPC codes must be in the following form: 



gi 




"10- 


•■ 


gim(t) 


g2 




1 ' 




g2m(t) 








•. 




gm—i 






•• 1 


gm—l,m (t) 


g m 




_ ••• • 


■ • 


gmm (t) 



where degg im < deg g mm = n — k. The information block u is represented as 
the vector u = (ui(t), ■ ■ ■ ,u m (t)), where 

" h-i 

J2 UijP i = 1, • • • ,m - 1, 

u m ,jt J i = m. 

j=n—k 

n—k—1 

The parity block u = (0, 0, ■ • • , 0, u m (t)), where u m (t) = Yl %n,j&, 1S the 

3=0 

remainder of u with respect to the reduced Grobner basis Q . The correspond- 
ing codeword is the result of subtracted vector u — u. This is received at the 
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Ul,0 ■ ■ ■ 



control signal 
Ui,0 ■ ■ ■ Ujjj-l 




r— | ( j— | 1 j storage of «i(t) out P ut 
►O^-t}^ HZ]"^" ! multiplying fflm(*) 



' 1 Uifi • • • Ui^-i 



storage of ltj(t) ► multiplying gi m (t) — 



^m,n—k ' ' ' ""m.lm— 1 



output 
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-o-i 



^"m,n—k ' ' ' 1^m,l m — l 



^m,0 ' ' ' ^777 , 71 — fe— 1 

output 



dividing by J mm (i) 



Figure 7: Serial-in serial-out architecture for type-II FG (EG and PG) LDPC 
codes. Input is information {uij}, and output is redundant parity bits {u m ,j}- 
The control signal is used to switch for feedback shift registers after entering 

Wl,0, • • • , Wl.Zi-l- 



output of architecture in Figure 0, serial-in serial-out architecture for FG 
LDPC codes. The element ® represents an adder (exclusive-OR element) 
and the rectangle represents a memory element (a shift register). The two 
remaining building elements correspond to multiplexer and gate elements. 
The gate element is a switch control with two status — open and close. The 
multiplexer element is signal choice control that selects signal either from 
input or from the feedback of shift registers. 

Then, the total number A m of adder elements for FG codes satisfies the 
following inequality: 

A m — Tn(n — k) < ml m < 2n. 
Moreover, the total number D m of required memory elements satisfies 

771—1 777 

D m < rnl m + y_ \k = (m - l)l m +}Jj< 2n. 

7=1 7 = 1 

Thus, we have proved that the hardware complexity of FG LDPC codes is 
0(n) order. 
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Table 1: Hardware complexity for several 3- dimensional type-II EG and PG 
LDPC codes. The first three rows evaluate type-II EG LDPC codes. The 
others evaluate type-11 PG LDPC codes. 
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n 


k 


n — k 


m 


adder 


memory 


1 


21 


15 


6 


3 


12 


26 


2 


315 


265 


50 


5 


76 


328 
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4599 


4227 


372 


9 


1681 


5769 


1 


35 


24 


11 


3 


16 


36 


2 


357 


296 


61 


5 


138 


438 


3 


4745 


4344 


401 


9 


1846 


6396 



For FG LDPC codes made from 3-dimensional EG and PG over the finite 
field F 2 s, where s = 1,2,3, we summarize computational results in Table 
[TJ The last two columns of Table [1] are the numbers of adder and memory 
elements, respectively. We see that the actual numbers of elements are less 
than the above estimation. 

7 Conclusions 

One contribution of this paper is to provide algorithms of computing Grobner 
basis for efficient systematic encoder of GQC codes. Our algorithms are appli- 
cable to not only binary GQC LDPC codes but also non-binary GQC LDPC 
codes and linear codes with nontrivial automorphism groups. Although the 
computation of Grobner basis is required only once at the construction of en- 
coder differently from decoding algorithm, our algorithms are still useful; for 
example, both algorithms can search effective codes rapidly in the polynomial 
(third power) order of code-length. For high-rate codes, we have shown that 
the algorithm applying transpose formula is faster than the echelon canonical 
form algorithm. It is expected that GQC LDPC codes improve the decoding 
performance of QC LDPC codes and make it close to that of the random 
LDPC codes. Another contribution of this paper is to demonstrate that the 
hardware complexity of the serial-in serial-out systematic encoder is the lin- 
ear order of code-length for FG codes and FG LDPC codes. By exploiting 
the structure of GQC codes, we believe that many new and optimum codes 
are constructed, and our results in systematic encoding might become a key 
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step to practical implementation. 
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